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Preface 


The  Photo-nuclear  Data  Center,  a  part  of  the  Center  for  Radiation  Research,  Is  one 
of  many  data  centers  within  the  National  Standard  Reference  Data  System  receiving  sup- 
port from  the  National  Bureau  of  Standards'  Office  of  Standard  Reference  Data,  Its 
function  is  to  collect,  abstract,  index,  and  evaluate  experimental  data  published  in 
the  world's  technical  literature  relating  to  the  interaction  of  electromagnetic  radia- 
tion with  nuclei.  It  provides  information-center  services  for  the  technological  and 
scientific  communities  concerned  with  photonuclear  reactions.  Ever  since  its  creation 
in  1963,  the  Data  Center  has  made  considerable  use  of  computer-compatible  index, 
bibliographic,  and  data  files.  These  are  used  with  the  central  computer  facilities  at 
the  National  Bureau  of  Standards  to  retrieve  specific  types  of  information;  to  prepare 
camera-ready  copy  for  publication  of  tabular,  bibliographical  and  graphical  information; 
and,  to  carry  out  various  manipulations  of  digital  data  sets  germane  to  their  evaluation. 

This  Technical  Note  is  the  first  of  several  being  written  to  describe  all  aspects 
of  the  computer-based  operations  of  the  Photonuclear  Data  Center.  These  notes  are  meant 
both  to  document  the  procedures  used  in  the  Data  Center  as  well  as  to  serve  as  'user's 
manuals'  for  individuals  having  a  need  to  make  use  of  the  Data  Center's  files.  The 
program  library  described  in  this  Note  has  been  designed  for  use  with  the  particular 
central  computer  facilities  available  at  the  National  Bureau  of  Standards.  No  attempt 
has  been  made  to  write  the  various  programs  in  a  completely  general  form  (facility- 
independent).  For  those  desiring  to  have  it,  the  complete  program  library  (in  symbolic 
form)  is  available  from  the  Photonuclear  Data  Center. 


Everett  G.  Fuller 
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The  NIRA  Computer  Program  Package 
(Photonuclear  Data  Center) 

H.  J.  Vander  Molen  and  Henry  M.  Gerstenberg 


The  Photonuclear  Data  Center's  NIRA  library  of  programs,  executable 
from  mass  storage  on  the  National  Bureau  of  Standards'  central  computer 
facility,  is  described.   Detailed  instructions  are  given  (with  examples), 
for  the  use  of  the  library  to  analyze,  evaluate,  synthesize,  and  produce 
for  publication,  camera-ready  tabular  and  graphical  presentations  of  digi- 
tal photonuclear  reaction  cross-section  data.  NIRA  is  the  acronym  for 
Nuclear  Information  Research  Associate. 


Introduction 


There  are  two  aspects  of  the  Photonu- 
clear Data  Center's  operations  which  make 
use  of  the  NBS  computer  facilities.  The 
first  is  related  to  the  abstracting  and  in- 
dexing of  experimental  papers  published  in 
the  world's  technical  literature.  This 
provides  the  input  both  for  the  Center's 
computer-searchable  index  to  the  available 
experimental  data  for  the  field  of  photonu- 
clear physics,  as  well  as  the  annotated 
bibliographies  or  data  indexes  which  the 
Center  regularly  publishes  and  updates  [l].1 

The  second  computer-oriented  activity 
of  the  Center  relates  to  its  library  of 
digital  cross  section  data  and  the  systema- 
tic evaluations  being  made  of  these  data. 
To  do  this,  it  has  been  necessary  to  gen- 
erate procedures  for  handling,  storing, 
classifying,  evaluating,  and  presenting 
data.  The  program  library  described  here 
is  an  implementation  of  these  objectives  on 
the  National  Bureau  of  Standards'  main  com- 
puter facility,  a  UNIVAC  1108  computer  op- 
erating under  the  EXEC-8  system.2 


2.  Description  of  the  Library 


The  NIRA  program  library  consists 
primarily  of  pr  Dgrams  which  handle  data, 


rather  than  index  or  bibliographic  informa- 
tion. Most  of  these  programs  are  quite 
simple,  both  conceptually  and  in  implementa- 
tion. All  can  be  executed  from  a  demand 
terminal,  although  ordinary  batch  usage  was 
the  primary  intention.  All  are  mutually 
compatible;  data  processed  by  any  program 
can  be  read  by  any  other  program  without 
human  manipulation.  It  is  likely  that,  in 
the  future,  need  will  arise  for  even  more 
programs.  However,  the  basic  operations 
described  here  should  suffice  for  most 
procedures . 

The  programs  reside,  as  addressable 
elements,  in  a  file  called  "NIRA*PROGM.". 
(The  asterisk  and  period  are  essential  parts 
of  the  name.)  This  file  is  stored  as  a 
"read-only"  file.  It  is  therefore  not  pos- 
sible to  alter  it  by  any  error  in  normal  use. 
(It  is  possible  to  delete  it,  however,  so  be 
cautious  when  using  DELETE  statements.)  Two 
or  more  runs  may  simultaneously  access  this 
file.  As  a  backup,  the  entire  program 
library  (in  symbolic  form),  as  well  as  a 
complete  listing,  is  kept  in  the  Photonucle- 
ar Data  Center.  At  present,  this  library 
deck  consists  of  more  than  3,000  punched 
cards  with  control  cards  in  place  for  re- 
loading to  mass  storage. 

Most  of  the  individual  programs  are 
designed  to  operate  on  a  deck  of  punched 
cards  which  define  a  "data  set".  Similarly, 
most  of  the  programs  output  either  real  or 


figures  in  brackets  indicate  the  literature  references  given  on  p.  6. 

2Use  of  brand  names  in  this  report  is  for  identification  purposes  only  and  is  not  to  be 
construed  as  a  recommendation  of  a  particular  product  or  company. 


virtual  punched  cards.  In  all  such  cases, 
the  program  adds  a  comment  card  to  the  deck 
which  states  what  the  program  has  done  to 
the  data.  This  card  is  listed  on  the  user's 
printout  for  the  run  In  which  the  program  Is 
called.  For  reasons  of  economy  (I.e.,  to 
limit  the  amount  of  the  computer's  memory 
used) ,  most  of  the  programs  will  accept  a 
maximum  of  500  data  points  and  100  comment 
cards.  This  can  be  changed  without  great 
difficulty,  but  this  quantity  of  memory  has 
been  more  than  sufficient  thus  far. 

As  defined  by  the  Photonuclear  Data 
Center,  a  "data  set"  consists  of  the  results 
of  the  measurement  of  an  excitation  func- 
tion, angular  distribution,  or  reaction 
product  spectrum  under  specified  conditions 
for  a  single  target.  The  set  consists  of  a 
Header  Card  which  contains  the  acquisition 
number,  number  of  data  points,  reference, 
nucleus,  and  reaction,  plus  some  data  source 
Information;  two  Label  Cards  which  give  the 
axis  labels,  usually  "Photon  Energy  (MeV)" 
and  "Cross  Section  (mb)";  an  arbitrary  num- 
ber of  Comment  Cards;  and  then  Data  Cards 
in  Fortran  format  (7(E9. V , ') ,E9.4,1X) . 


A  few  programs  require  more  elaborate 
control  cards  than  this  MAP-IN-XQT  sequence. 
These  have  been  combined  with  the  appropri- 
ate individual  programs  into  what  in  the 
following  discussion  are  referred  to  as 
"job-control  decks".  Each  deck  is  named, 
and  only  one  card  need  be  entered  to  call  a 
deck;  e.g.,  PPLOT  is  called  by: 

@ADD,P  PF. PPLOT 

The  individual  writeups  give  the  exact 
cards  to  use.   (Do  not  forget  the  ASG-USE 
sequence  at  the  beginning  of  your  run, 
however ! )  Note  that  the  run  cards ,  tape 
identification  labels,  etc.,  that  are  used 
in  a  number  of  the  examples,  are  given  for 
purposes  of  illustration  only,  and  that  In 
making  use  of  the  library,  the  user  should 
be  sure  to  use  the  appropriate  number  and 
designation  associated  with  his  computer 
account.  Throughout  the  listings,  square 
brackets,  i.e.,  [  ],  are  used  to  designate 
cards  containing  information  germane  to 
particular  runs  only.  Brackets  are  not 
used  in  entering  information  unless  speci- 
fically called  for  in  the  program  descrip- 
tion. 


3.  Use  of  the  Program  Library 


To  use  the  library,  it  is  necessary  to 
assign  the  library  file  to  a  particular 
computer  name.  To  do  this,  use  the  two 
EXEC~B  control  cards  [2]: 

@ASG,A  NIRA*PROGM. 

@USE  PF.,NIRA*PROGM. 

The  "commercial  at"  sign  is  represented  by 
a  7/8  punch  on  the  card  punch,  and  by  "@" 

on  the  Teletype. 

As  indicated  in  the  individual  descrip- 
tions given  in  the  Appendix,  most  of  the 
programs  are  called  by  entering  three  cards 
For  example:  SORT  is  called  by: 

@MAP,INX  MAPR,S0RT 

IN  PF.SORT 
@XQT  SORT 

Similarly,  SMOOTH  is  called  by: 

@MAP,INX  MAPR, SMOOTH 

IN  PF.  SMOOTH 
@XQT  SMOOTH 


Qualitative  Description  of  the  Programs 
and  Their  Use 


A  brief  description  of  the  purpose  of 
the  various  programs  is  given  here.  The 
individual  descriptions  should  be  consulted 
before  they  are  used. 


4.1.  Programs  which  Create  Data  Sets 

DIGIT  is  a  job-control  deck  that  takes 
data  from  a  tape  made  by  the  sonic  delay- 
line  digitizer,  converts  it  to  data-library 
format,  operates  on  it  with  CLEAN,  provides 
punched  cards  for  the  data  library,  and/or 
leaves  the  data  in  a  computer  file  element 
called  DIGOUT  for  use  by  other  programs. 

HIST  converts  a  normal  data  set  into  a 
histogram. 

UNHIST  converts  a  histogram  into  a 
normal  data  set . 


4.2.  Programs  Which  "Dress-Up"  Data 

SMOOTH  performs  a  normal  vertical 
smoothing  on  raw  data. 

STIFF  is  used  to  "dress-up"  digital 
data  so  that  it  may  be  used  as  input  to  ob- 
tain computer-drawn,  smooth,  continuous 
cross  section  curves. 

AVRG  smooths  raw  data  over  a  pre-set 
energy  interval. 

INTERP  parabolically  interpolates 
smoothed  data  to  give  points  at  "nice" 
energies . 

SQUARE  squares  up  histograms. 


4.3.  Programs  Which  Manipulate  Data 

COMBIN  algebraically  adds  two  cross 
sections.  It  will  read  a  (y,n)  and  a 
(Yj2n)  cross  section,  and  punch  a  (y>xn) 
cross  section,  for  example. 

DBAL  reverses  a  cross  section  by  means 
of  the  principle  of  detailed  balance. 

CALC  renormalizes  data  sets. 

UNSTAT  removes  the  statistical  theory's 
multiplicity  correction  from  a  (Y5sn)  cross 
section. 


4.4.  "Housekeeping"  Programs 

CLEAN  checks  data  for  format  and  makes 
several  computer-oriented  (non-manipulative) 
corrections.  A  data  set  which  has  passed 
through  CLEAN  should  give  no  problems  when 
read  by  the  other  programs. 

MONO  forces  energies  to  increase 
monotonically . 

SORT  sorts  data  points  by  energy. 

TRUNC  picks  out  a  subset  of  a  data  set. 

CONCAT  joins  together  several  data 
sets  covering  different  energy  ranges. 

4.5.  Programs  Which  Produce  Final  Copy 

YTELD4  is  used  to  print  a  column-wise 
table  of  cross  sections.  It  is  also  used 
to  calculate  c_,  at  10,  15,  20,  25,  and 
30  MeV.       x 


SPCDST  is  a  deck  used  to  print  an 
index  table  of  spectral  and  angular  distri- 
bution data. 

LIST  is  used  to  produce  camera-ready 
copy  of  tables  made  by  either  YIELD4  or 
SPCDST. 

PPLOT  is  a  plotting  job-control  deck. 
It  produces  plots  which  are  suitable  for 
mounting  on  reduction  paper. 

GDR  calculates  the  integral  moments  a0, 
a_i,  and  o_2  of  the  cross  section,  and  esti- 
mates the  giant  resonance  parameters  E  ,  a  , 
and  r.  See  Ref.  [1]  for  definitions.  m   m 

4.6.  Use  of  Programs 

The  use  of  the  programs  is  completely 
flexible.  Nevertheless,  there  is  a  "normal" 
or  "usual"  sequence  of  steps  through  which 
data  pass  from  journal  article  to  final 
copy. 

(1)  The  data  are  digitized  and 
assigned  an  acquisition  number.  The  data 
set  is  then  run  through  DIGIT.  If  the  data 
already  exist  in  digital  form,  they  are 
punched  and  run  through  CLEAN.  UNHIST  may 
also  be  used  if  the  data  are  given  as  a 
histogram. 

(2)  If  CLEAN  indicates  monotonicity 
problems  with  the  data,  SORT  and/or  MONO 
may  be  used  to  correct  the  trouble. 

(3)  The  data  set  is  manipulated,  if 
appropriate,  to  get  it  into  the  proper  units 
and  form  for  comparison.  COMBIN,  DBAL, 
UNSTAT,  TRUNC,  and  CALC  are  provided  for 
such  manipulation. 

(4)  The  data  are  plotted,  using  the 
CPLOT  mode  of  PPLOT.  Comparisons  are  made, 
and  the  best  data  are  selected. 

(5)  The  selected  data,  if  a  (y.sn), 
(y,xn),  or  (y,y,  ),  are  run  through  GDR  to 
calculate  a„ ,  a  '  ,  and  a_„ . 

also  gives  estimated  values  for  E 

n 
and  T. 


This  program 


(6)  The  selected  data  are  used  to 
create  a  smoothed  cross  section  curve.  This 
smooth  curve  is  the  source  of  the  evaluated 
data  set.  The  curve  is  usually  drawn  by 
hand.  SMOOTH  and  AVRG  are  sometimes  used 
in  conjunction  with,  or  even  in  place  of, 
the  manual  drawing.  (SQUARE  is  used  to 
dress  up  histograms . )  If  the  manual  method 
is  used,  the  hand-drawn  line  is  digitized, 


run  through  DIGIT  and  MONO,  and  then  run 
through  STIFF  to  smooth  out  the  hand-induced 
irregularities.  If  the  hand-drawn  line  was 
digitized  in  sections,  CONCAT  is  used  to 
join  the  sections  together  before  STIFF  is 
used. 

(7)  The  evaluated  data  are  then  run 
through  INTERP  to  obtain  a  compact  repre- 
sentation with  cross-section  values  given  at 
specified  energy  intervals.  Occasionally, 
MONO  is  needed  between  STIFF  and  INTERP. 
The  output  of  INTERP  is  filed  as  the  final 
form  of  the  "evaluated"  cross  section. 

(8)  The  yields  are  calculated  from  the 
selected  (not  evaluated)  data  using  the 
YIELD  aspects  of  YIELD4. 

(9)  The  evaluated  cross  sections,  as 
well  as  the  selected  data  set,  are  plotted 
using  PPLOT. 

(10)  The   evaluated  cross  sections  are 
tabulated  using  YIELDS  and  LIST." 

(11)  The  index  to  supplementary  data 
is  tabulated  using  SPCDST. 


5.  Some  Notes  on  Efficient  Use 


Almost  all  of  the  programs  described 
in  the  preceding  pages  read  cards  for  input 
and  punch  cards  for  output.  It  is  not  nec- 
essary to  execute  just  one  program  per 
computer  run.  There  are  several  "tricks" 
which  enable  the  user  to  perform  a  sizable 
amount  of  work  in  one  computer  run. 


5.1.  Multiple  Execution 

Most  of  the  programs  are  called  in  the 
following  manner: 

@MAP,INX  MAPR,GDR 

IN  PF.GDR 
@XQT  GDR 

The  first  two  cards  direct  the  MAP 
processor  to  transform  the  relocatable 
element  PF.GDR  into  the  absolute  (execut- 
able) element  GDR.  Even  after  execution, 
the  absolute  element  GDR  is  still  in  the 
run's  temporary  program  file  (TPF$.).  It 
can  be  re-executed  with  another  XQT  card;. 


@MAP,INX  MAPR,GDR 

IN  PF.GDR 
@XQT  GDR 

[1st  data  set]  i 

@XQT  GDR 
[2nd  data  set] 
@QXT  GDR 
[3rd  data  set] 
etc. 

5.2.  Temporary  File 

A  temporary  file  is  useful  when  the 
same  data  set  is  needed  for  input  to  more 
than  one  program.  The  file  is  created  by: 

@ASG,T  TF. 

Data  are  then  put  into  the  file  element  by 
element .  The  file  elements  are  fed  to  the 
processing  programs  by  ADD. 

Example:  Suppose  one  started  out  with  a 
(y3n+np)  and  a  (y,2n)  data  set.  .  It  is 
desired  to  calculate  a  (y^sn)  and  a  (Y,xn)' 
data  set.  The  cards  would  be: 

@ASG,T  TF. 

<§ELT,IS  TF.l 

[(y,n+np)  data]  ., 

@ELT,IS  TF.2 

[(y,2n)  data] 

@MAP,INX  MAPR,COMBIN 

IN  PF.COMBIN 
@XQT  COMBIN 
0,0,1. ,1., 
@ADD,P  TF.l 
@ADD,P  TF.2 
@XQT  COMBIN 
0,0,1. ,2., 
@ADD,P  TF.l 
@ADD,P  TF.2 


5.3-  Breakpoint  File 

A  breakpoint  file  is  useful  when  it  is 
desired  to  feed  the  output  of  one  program 
into  the  input  of  another  program.  The  file 
is  created  by: 

@ASG,T  BP. 

The  punched  output  of  a  program  is  diverted 
to  the  breakpoint  file  by: 

@BRKPT  PUNCH$/BF  (Note  no  punctuation!) 

Then  the  first  program  is  executed.  After 
the  last  data  card  of  the  first  program,  the 
punch  channel  is  re-activated  by: 

SBRKPT  PUNCH$ 

Now  the  contents  of  the  breakpoint  file  can 
be  fed  into  the  2nd  program  by  means  of: 

@ADD,P  BF. 

Example:  Suppose  it  is  desired  to  calcu- 
late  a  (y,xn)  total  neutron  yield.  The 
starting  data  are  a  (Y>sn)  data  set 
measured  via  a  BF3  detector  and  a  statisti- 
cal multiplicity  correction.  It  is  neces- 
sary to  first  call  UNSTAT  to  transform  the 
(YjSn)  data  set  into  a  (y,xn)   data  set, 
then  to  call  YIELDM  to  calculate  the  yield. 


@ASG,T  BF. 

@MAP,INX  MAPR,  UNSTAT 

IN  PF. UNSTAT 
ffiRKPT  PUNCH$/BF 
@XQT  UNSTAT 
[parameter  card] 
[(YjSn)  data] 
@BRKPT  PUNCH$ 
@MAP,INX  MAPR, YIELDS 

IN  PF.YIELD4 
@XQT  YIELD4 

(NONE)  YIELD 
@ADD,P  BF. 

The  breakpoint  file  is  erased  whenever 
new  gBRKPT  PUNCH$/BF  card  is  encountered, 
iesides  adding  data  with  an  @ADD  statement, 
he  contents  of  the  breakpoint  file  can  be 
isted  by  entering: 

@ADD,P  PF.LIST 


as  can  be  seen  in  the  description  of  SPCDST. 
The  breakpoint  file  can  also  be  punched  by: 

@ELT,ID  PUNCH 

6ADD,D  BF. 

@END 

@PCH,S  PUNCH 

(The  SYM  processor  of  EXEC-8  will  also  punch 
BF,  but  requires  cataloging.  This  had  led 
to  interference  between  runs . ) 

The  breakpoint  file  is  an  extremely 
powerful  tool  when  used  in  conjunction  with 
a  temporary  file.  The  following  example 
begins  with  a  (p,Y0)  data  set.  The  data  set 
is  inverted  and   plotted.  Simultaneously, 
a  smoothed  form  of  the  same  cross  section  is 
digitized,  smoothed,  punched,  inverted, 
plotted,  and  tabulated. 


@RUN,M/NR  VAND41,30850-VANDER,NIRA, 1,50/1000 

@ASG,A  NIRA*PROGM. 

@USE  PF.,NIRA*PROGM. 

@ASG,T  BF. 

@ASG,T  TF. 

@MAP,INX  MAPR,DBAL 

IN  PF.DBAL 
gBRKPT  PUNCH$/BF 
@XQT  DBAL 
[parameter  card] 
[(p,YQ)  data] 
(3BRKPT  PUNCH$ 
@ELT,IS  TF.l 
@ADD,P  BF. 
@ELT,IS  DIGIN 
[digitizer  control  cards] 
@ASG,LJ0  IN.,8C,VAN02 
@SETC  1 

@ADD,P  PF. DIGIT 
@MAP,INX  MAPR, MONO 

IN  PF.M3NO 
@BRKPT  PUNCH$/BF 
@XQT  MONO 
@ADD,P  DIGOUT 
@BRKPT  PUNCH$ 
@ELT,IS  TF.2 


@ADD,P  BF. 
@MAP,INX  MAPR, STIFF 

IN  PF. STIFF 
gBRKPT  PUNCH$/BF 
@XQT  STIFF 
[parameter  card] 
@ADD,P  TF.2 
@BRKPT  PUNCH$ 
@ELT,IS  TF.2 
@ADD,P  BF. 
@BRKPT  PUNCH$/BF 
@XQT  DBAL 
[parameter  card] 
@ADD,P  TF.2 
@BRKPT  PUNCH$ 
@ELT,IS  TF.2 
@ADD,P  BF. 
@MAP,INX  MAPR,INTERP 

IN  PF.INTERP 
(9BRKPT  PUNCH$/BF 
@XQT  INTERP 
@ADD,P  TF.2 
@BRKPT  PUNCH$ 


@ELT,IS  TF.2 
@ADD,P  BF. 
@PCH,S  TF.2 
@MAP,INX  MAPR,YIELD4 

IN  PF.YIELD4 
(2BRKPT  PUNCH$/BF 
@XQT  YIELD4 
[blank  card] 
@ADD,P  TF.2 
gBRKPT  PUNCH$ 
@ADD,P  PF.LIST 
@ADD,P  PF.PPLOT 
@ADD,P  TF.l 
$CURVE 
@ADD,P  TF.2 
$PPLOT 

@ADD,P  PF.PRINTR 

@MSG,W  30850-VANDER,WRITE  ENABLE  REEL  VAN03 
@ASG,TMI  9 . ,8C,VAN03W 
@CLOSE  9. 
@XQT  NBS*CALCMP. 

@FIN 
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AVRG 


6.  Appendix;  NIRA  Listings 
6.1.  Programs 
AVRG  is  a  means  of  smoothing  a  data  set.  The  cross  section  at  each  energy,  E.  is 


replaced  by  the  average  value  over  a  specified  energy  interval,  e.  (This  is  in 

contrast  to  SMOOTH  which  averages  over  3  points  regardless  of  spacing.)  The  averaged  data 

set  is  then  punched.  The  formula  used  is: 

i       A+  £/2 

a(E.)  =  i    /  x  a(E)dE  . 


l    e 


r 

J  E.- 


E.-  e/2 


The  integration  is  rectangular.  The  averaging  interval  "e"  is  constant  except  at  the 
low  energy  end.  If  E.-  E  <  e/2,  e  =  2(E.-  E  ). 

The  averaging  interval  "e"  is  punched  as  a  floating  point  number  in  free  form, 
followed  by  a  comma,  on  a  parameter  card  immediately  preceding  the  data  set.  If  the  card 
is  omitted,  e  =  1  MeV.  Any  number  of  data  sets  may  be  entered.  The  value  given  to 
epsilon  remains  in  force  until  a  new  parameter  card  is  encountered. 

@MAP,INX  MAPR,AVRG 

IN  PF.AVRG 
@XQT  AVRG 
[parameter  card] 
[data  set] 

CALC 

CALC  manipulates  a  data  set  by  the  formulas:  E.=  A  E.+  A    and  a.=  B  a  +  B„  .  The 

first  card  read  by  CALC  must  contain  A,,  A„,  B,  ,  B?,  and  the  new  acquisition  number.  If 

zero  is  entered  for  the  new  acquisition  number,  the  old  acquisition  number  will  be  used. 
The  program  will  read  only  one  data  set . 

Example:  If  it  is  desired  to  change  a  data  set  scale  from  millibarns 

to  microbarns,  the  parameter  card  would  be:  1. ,0. ,1000. ,0. ,0 

@MAP,INX  MAPR,CALC 

IN  PP. CALC 
@XQT  CALC 
[parameter  card] 
[data  set] 

CLEAN 

CLEAN  reads  a  data  set  and  checks  it  for  format.  It  will  automatically  modify  the 
set  as  follows: 


(1)  If  the  label  cards  are  missing,  "PHOTON  ENERGY  (MEV)"  and  "CROSS  SECTION  (MB)" 
are  inserted. 

(2)  If  the  number  of  data  points  is  missing  or  incorrect,  the  correct  number  is 
inserted. 

(3)  If  the  header  card  is  missing,  one  is  inserted. 

(4)  If  two  adjacent  data  points  have  the  same  energy,  the  second  energy  is  increased 
one  part  per  thousand. 

(5)  If  the  data  are  non-monotonlc,  except  for  the  case  above,  a  warning  is  printed. 

(6)  If  a  blank  card  is  encountered,  it  is  skipped  over. 

(7)  If  two  points  have  identical  x-  and  y-coordlnates  (e.g.,  due  to  double-sparking 
of  digitizer),  the  second  point  of  the  doublet  is  deleted. 

This  routine  can  read  any  "free"  format.  A  warning  Is  printed  if  the  data  are  not  in  the 
standard  "(7(E9.4,', ' ),E9.4,1X)"  Data  Center  format.  The  monotonicity  and  double-sparking 
tests  are  bypassed  if  the  acquisition  number  is  negative  (i.e.,  special  test  data). 

A  data  set  punched  by  CLEAN  can  be  handled  by  all  the  other  programs  without  problems. 
@MAP,INX  MAPR,  CLEAN 

IN  PF. CLEAN 
@XQT  CLEAN 

COMBIN 

COMBIN  is  used  to  algebraically  add  two  cross  sections  and  punch  a  new  data  set.  The 
formula  is: 

c'  =  aa  +  6a,  , 
a    b  ' 

where  a  and  8  can  be  any  positive  or  negative  floating  point  number.  The  two  data  sets 
must  be  preceded  by  a  parameter  card  containing  NA,  NB,  a,  and  8  in  free  field  format.  NA 
and  NB  must  be  fixed-point  positive  or  negative  integers,  and  a  and  8  are  floating-point 
numbers.  The  maximum  energy  in  the  new  data  set  will  be  the  larger  of  that  in  the  two  sets 
combined  if  NA  =  0.  Setting  NA  =  -1  results  in  the  maximum  energy  being  the  lesser  of  that 
in  the  two  sets.  In  both  cases,  NB  is  the  acquisition  number  for  the  new  data  set. 

Only  two  data  sets  may  be  entered.  The  set  with  the  best  resolution  should  be  the 
second  data  set.  If  the  resolution  differs  greatly  between  the  two  sets,  then  the  programs 
AVRG  or  SMDOTH  should  be  used  to  obtain  comparable  resolution  before  using  COMBIN. 

The  program  makes  the  resolution  of  output  data  similar  to  that  of  the  first  data  set. 
For  the  region  where  the  data  sets  overlap  in  energy,  the  second  set  is  linearly  interpo- 
lated to  match  exactly  the  energies  of  the  first  data  set.  For  the  regions  where  there  is 
no  overlap  and  data  exist  only  in  the  second  set,  then  the  resolution  of  the  first  set  is 
averaged  over  ten  data  points  nearest  the  region  of  interest.  This  average  energy 
resolution  is  then  used  in  the  linear  interpolation  of  the  second  data  set  to  obtain  output 
data  points  of  the  second  data  set  having  a  resolution  similar  to  that  of  the  first  data 
set. 

@MAP ,  LNX  MAPR,  COMBIN 

IN  PF.  COMBIN 
@XQT  COMBIN 
[parameter  card] 
[data  set  A] 
[data  set  B] 


CONCAT 

This  program  reads  an  arbitrary  number  of  data  sets,  and  punches  one  data  set  which 
is  a  simple  concatenation  (i.e.,  joining  together)  of  the  input  data  sets.  The  header  and 
label  cards  of  the  output  are,  except  for  the  number  of  points,  copied  from  those  of  the 
first  data  set  read.  The  total  number  of  points  must  not  exceed  500.  Generally,  it  is 
wise  to  feed  the  output  of  CONCAT  into  SORT. 

@MAP , INX  MAPR, CONCAT 

IN  PF. CONCAT 
@XQT  CONCAT 
[data  sets] 


DBAL 

DBAL  reverses  a  cross  section  via  the  principle  of  detailed  balance.  It  reads  a  para- 
meter card  and  a  data  set. 

For  reactions  of  the  type  A+y  rf  B+b,  the  equations  used  are: 


a(Y,bQ)  =  f(A,Ee)  o(b,Y0)     where, 


The  parameters  are: 


(2J  +l)(2s+l)  m  Mg  (E  -BE) 

f(A,E  )  =  7 „ T  -  =  t ..  ,   with 

e       (2JA+1)     Mg+m    £2 


E  =  =5-2-  E,+  BE   or   E  =  E  -  * — 

e    Mg+m   b  e    Y    2(M  2+E  ) 

A    y 


q  =  flag  to  indicate  direction  of  change 
m  =  mass  of  particle  b 
s  =  spin  of  particle  b 
M.  =  mass  of  nucleus  A 


% 


=   mass  of  nucleus  B 


Jfl  =  spin  of  nucleus  A 

J„  =  spin  of  nucleus  B 
ts 

BE  =  binding  energy  of  b  in  A 

To  turn  a  ground-state  capture  cross  section  into  the  related  photodisintegration 
cross  section,  the  parameter  card  is: 


1  ,  m,  s,  Mg,  JB,  MA,  JA, 


BE 


Similarly,  for  the  reverse  transformation,  the  card  should  read: 

0.,  m,  s,  Mg,  JB,  MA,  JA,  BE 

Any  number  of  data  sets  may  be  entered,  but  each  must  have  a  parameter  card. 

@MAP,INX  MAPR,DBAL 

IN  PF.DBAL 
@XQT  DBAL 
[parameter  card] 
[data  set] 

GDR 

GDR  reads  a  data  set,  calculates  an,  a   ,,  and  a   „,  and  estimates  E_,  a   ,  and  the 

0-1      -2'  0'  max' 

width.  Only  one  data  set  may  be  entered. 

The  integration  is  done  by  the  trapezoidal  rule.  The  upper  limit  of  integration  is 
the  greatest  integer  which  is  less  than  the  max  energy  of  the  data  set.  The  value  of  E_ 

is  a  calculated  mean  energy.  The  estimates  of  o    and  width  are,  in  practice,  almost 
worthless.  A  better  job  can  be  done  by  eye. 

The  data  set  is  repunched  with  the  calculated  giant  resonance  parameters  given  on  a 
comment  card. 

@MAP,INX     MPR,GDR 

IN     PF.GDR 
@XQT     GDR 
[data  set] 

HIST 

This  program  converts  a  normal  data  set  into  a  histogram,  vertical  bars  are  placed 
midway  between  every  two  points  of  the  input  data  set.  A  slight  inclination  of  the 
verticals  is  Introduced  to  preserve  monotonicity.  Any  number  of  data  sets  may  be  entered. 
CAUTION  -  The  output  may  look  strange  if  the  input  is  not  approximately  evenly  spaced  in 
energy . 

@MAP,INX     MAPR,RTST 

IN  PF.HIST 
@XQT  HIST 
[data  sets] 

INTERP 

INTERP  reads  a  data  set  and  punches  a  parabolic  interpolation  of  the  input .  The 
output,  except  for  the  first  and  last  points,  is  given  at  every  tenth,  fifth,  or  half  MeV, 
depending  on  local  curvature. 

The  output  is  suitable  for  input  into  PPIOT  and  YIELDS.  Any  number  of  data  sets  may 
be  entered. 
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@MAP,INX  MAPR,INTERP 

IN  PF.INTERP 
@XQT  INTERP 
[data  sets] 


MONO 


1VDN0  is  a  program  that  corrects  monotonicity  problems  which  result  from  roundoff  er- 
rors in  going  from  the  raw  digitized  data  to  the  Data  Center's  data  card  format.  The  order 
of  the  points  is  maintained,  but  energies  are  adjusted  upwards  such  that  they  monotonically 
increase.  (This  is  not  always  the  appropriate  correction  -  see  SORT.) 

If  more  than  one  data  set  is  entered,  each  one  will  be  checked,  corrected  if  neces- 
sary, and  punched. 

@MAP,INX  MAPR, MONO 

IN  PF.MONO 
@XQT  MONO 
[data  sets] 

SMOOTH 

SMOOTH  performs  one  pass  of  a  normal  vertical  smoothing  on  a  data  set.  Any  number  of 
data  sets  may  be  entered.  The  energies  do  not  need  to  be  equally-spaced.  The  algorithm 
used  is: 


1 
°i  =  2 


a     a. 

a.  +  _1+1  ~  i-J-  (e.  -  E.  ,)  +  a.  , 

1    E.  ,,  -  E.  ,   1    l-l     l-l 
l+l    l-l 


In  the  special  case  of  equally-spaced  energies,  this  reduces  to  the  familiar  form: 

1      1.1    .  1 

o.  =7-a.  -,  +  -a.  +T-O.J, 
i   4  i-I   2  i   4  i+I 

@MAP , INK  MAPR, SMOOTH 

IN  PF. SMOOTH 
@XQT  SMOOTH 
[data  sets] 

SORT 

This  program  sorts  the  points  of  a  data  set  by  energy.  Its  primary  uses  are: 
(1)  removing  monotonicity  problems  caused  by  digitizing  data  in  a  region  of  high-point 
count  density,  and  (2)  fixing  a  deck  of  cards  someone  dropped  on  the  floor.  Note  that 
SORT  and  MONO  have  the  same  goal  but  accomplish  it  in  different  ways;  SORT  changes  the 
order  of  the  points,  while  MONO  changes  the  value  of  the  points. 

If  more  than  one  data  set  is  entered,  the  program  will  sort  each  data  set. 

@MAP,INX  MAPR, SORT 

IN  PF.SORT 
@XQT  SORT 
[data  sets] 
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SQUARE 

This  program  is  used  to  "dress  up"  a  histogram-type  bar  graph.  (We  use  "histogram" 
to  denote  a  bar  graph,  not  in  the  strict  statistical  sense.)  If  the  "corner  points"  of  a 
histogram  are  digitized: 


~L 

T — T 

• • 

f — • 

n '> 

then,  in  general,  the  output  of  the  digitizer  program  will  give  a  graph  in  which  the  hori- 
zontal and  vertical  lines  are  not  truly  parallel  to  the  axes. 

SQUARE  reads  such  a  "raw"  data  set,  and  adjusts  the  points  such  that  horizontal  lines 
are  truly  horizontal.  Vertical  lines  are  inclined  very  slightly  (one  part  per  thousand) 
to  avoid  monotonicity  problems. 

The  height  of  each  "squared"  horizontal  is  the  average  of  the  heights  of  the  two  raw 
points  defining  it,  and  similarly  for  the  verticals.  Any  number  of  data  sets  may  be 
entered . 

@MAP , INX  MAPR, SQUARE 

IN  PF. SQUARE 
@XQT  SQUARE 
[data  sets] 

STIFF 

STIFF  performs  one  pass  of  a  non-vertical  smoothing  routine.  This  means  that  energies 
are  also  adjusted  to  smooth  the  data.  Such  smoothing  is  the  appropriate  treatment  for 
digitized  curves  where  errors  are  both  vertical  and  horizontal. 

The  program  is  used  to  smooth  digital  data  which  has  been  generated  from  "smooth" 
curves  drawn  by  hand  through  actual  data  points.  The  resulting  "smooth"  data  are  used  to 
generate  computer-drawn  curves  which  approximate  the  actual  data  presented  in  graphical 
form.  They  are  also  used  to  generate  the  tables  of  "smoothed"  cross  section  data. 

The  smoothing  is  carried  out  by  averaging  over  three  successive  data  points,  i,  j,  and 
k,  when  their  coordinates  are  expressed  in  inches.  The  program  first  uses  points  i  and  k 
to  determine  a  straight  line.  The  coordinate  axes  are  then  rotated  to  establish  an  x'-axis 
which  is  parallel  to  this  line,  and  the  average  y' -coordinate  of  the  point  j  is  determined 
from: 

y'(j)  =  Jy'(i)  +|y'(J)  +  Jy'OO. 
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The  coordinates  of  the  point  j,  x'(j),  y'(j)  are  then  referred  back  to  the  original  un- 
rotated  system  to  obtain  the  "smoothed"  coordinates  of  the  point  j. 

Since  the  program  must  work  with  the  coordinates  of  each  point  expressed  in  actual 
linear  dimensions ,  the  first  card  must  contain  XLNGTH  and  YLNGTH  ( in  inches ) ,  and  the  three 
fiduciary  points  of  the  plot  from  which  the  data  were  digitized. 

Example: 


-i-  80 


60 


in  40 

in 


20 


i    o 


20 


30 


6.0- 


40 


The  parameter  card  is  then: 


6.,5.5,10.,0.,40.,0.,10.,* 


More  than  one  data  set  may  be  entered.  If  new  parameter  cards  are  not  entered  along 
with  the  new  data  sets,  the  old  parameters  will  remain  in  effect. 

The  primary  use  of  STIFF  is  to  "dress  up"  a  smooth  curve  by  reducing  the  small  point- 
to-point  irregularities  introduced  by  hand-digitizing.  The  output  is  suitable  for  input 
into  INTERP. 

@MAP,INX  MAPR, STIFF 
IN  PF. STIFF 
@XQT  STIFF 
[parameter  card] 
[data  sets] 
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TRUNC 

This  program  reads  a  data  set  and  punches  a  truncated  data  set  consisting  of  data 
points  which  fall  within  a  selected  energy  interval.  The  upper  and  lower  limits  are 
entered  on  a  parameter  card  as  floating  point  numbers  in  free  format.  The  energy  interval 
is  inclusive  of  the  end  points.  Any  number  of  data  sets  may  be  entered.  If  no  new  energy 
limits  are  entered,  the  old  limits  remain  in  force.  If  no  limits  are  entered  [0.,30.,] 
is  assumed. 

@MAP,INX  MAPR, TRUNC 

IN  PF.TRUNC 

@XQT  TRUNC 

[E,    ,  E     ] 
T_ower'  upper 

[data  sets] 

UNHIST 

This  program  reads  a  histogram  (either  "squared"  or  "raw")  and  converts  it  into  a 
normal  data  set. 

Output  points  are  placed  in  the  center  of  the  horizontal  bars,  except  when  multiple 
bins  are  detected.  No  problem  should  occur  if  (as  is  almost  always  the  case)  the  input 
bins  are  of  approximately  equal  width.  Any  number  of  data  sets  may  be  entered. 

@MAP,INX  MAPR, UNHIST 

IN  PF. UNHIST 
@XQT  UNHIST 
[data  sets] 

UNSTAT 

UNSTAT  reads  a  (y5sn)  cross  section,  and  calculates  a  (y>xn)  cross  section  via  the 
statistical  model.  Its  primary  use  is  to  strip  out  the  statistical  multiplicity  correction 
used  by  many  experimenters  with  Halpern  detectors.  The  relations  used,  taken  from 
T.  K.  Deague,  E.  G.  Muirhead,  and  B.  M.  Spicer,  Nucl.  Phys.  A139,  501  (1969),  are: 

a(sn)  =  a(y5n)  +  a(y3np)  +  a(y52n)  +  

a(xn)  =   a(y,n)  +  a(y,np)  +  2a(y,2n)  +  

a(sn)  =  a(xn)  F(E) 

F(E)  =  1+  (l-x)Cl-  (l-f(E))exp-f(E)] 


f(E)  =  (E-E2)  (a/(E-El))1/2,       where 


x  =  direct  de-excitation  fraction 
a  =  level  density  (in  MsV-1) 
El  =  the  (y,n)  threshold  (in  MeV) 
E2  =  the  (y,2n)  threshold  (in  NfeV) 
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The  first  card  must  contain  floating-point  numbers  for  x,  a.  El.  and  E2  in  free  format. 

Any  number  of  data  sets  may  be  entered.  If  new  parameter  cards  are  not  entered  along 
with  the  new  data  sets,  the  old  parameters  will  remain  in  effect. 

@MAP,INX  MAPR,UNSTAT 

IN  PF.UNSTAT 
@XQT  UNSTAT 
[parameter  card] 
[data  sets] 

YIELDS 

The  program  YIELDS  is  designed  to  operate  on  a  number  of  data  sets  and  calculate: 
(1)  a  cross-section  table  with  values  for  each  of  the  reactions  entered  being  given  at 
evenly-spaced  energy  intervals  (usually,  every  half  MeV) ,  and/or  (2)  values  of  a_, 
for  upper-energy  integration  limits  of  10,  15,  20,  25,  and  30  MeV  for  each  of  the 
reactions  entered.  A  linear  interpolation  is  used  to  calculate  values  at  specified  energy 
values.  This  program  is  a  version  of  YIELD3  which  operates  on  a  single  data  set  to 
produce  a  printout  in  which  for  each  energy  E.  in  the  set  the  following  quantities  are 
given : 


E. 

l 


a./E. 


a./E? 
l  l 


i 


X>./e. 


p/Ej 


The  operations  carried  out  by  YIELDS  are  determined  by  entries  punched  on  a  parameter 
card.  These  entries  must  be  separated  by  at  least  one,  and  no  more  than  five,  blank  spaces. 
There  must  be  no  blanks  within  an  entry,  and  commas  must  not  be  used  to  separate  entries. 
This  parameter  card  must  be  provided  even  if  no  entries  are  punched  on  it.  In  its  normal 
state,  the  program  calculates  a  cross-section  table  at  half-MeV  intervals.  The  following 
entries  may  be  made: 


DELTA=XXXXXX 


(G,XN) 


62Ful 


PUNCH 


(NONE) 


If  an  interval  other  than  0.5  MeV  is  desired, 
replace  the  six  X*s  with  a  floating-point 
number  giving  the  interval.  When  the  normal 
DELTA  is  not  used,  this  must  be  the  first  par- 
ameter on  the  card. 

To  skip  a  particular  reaction  in  building  the 
table,  enter  the  reaction  with  parentheses. 

To  skip  a  particular  reference  in  building  the 
cross-section  table,  enter  the  reference  number. 

The  output  is  listed  in  the  normal  Data  Center 
cross-section  library  format,  i.e.,  Energy, 

Cross  Section,  Energy,  Cross  Section etc. 

Note:   "PUNCH"  and  "YIELD"  are  mutually  exclusive. 

No  cross-section  table  is  desired. 
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YIELD  The  yield  table  is  calculated,  i.e.,  c_,  values 

at  10,  15,  20,  25,  and  30  MeV. 
Note:  "PUNCH"  and  "YIELD"  are  mutually  exclusive. 

PRINT  Gives  a  full  printout  of  YIELD3  for  each  data  set 

entered.  See  above  discussion. 

The  cross-section  table  is  punched  in  the  proper  format  to  be  listed  by  LIST.  See  the 
description  of  LIST  for  more  details.  Note:  The  cross-section  table  can  hold  a  maximum 
of  200  energies  and  13  reactions. 

@MAP,INX  MAPR, YIELDS 

IN  PF. YIELDS 
@XQT  YIELDS 

[parameter  card]     This  card  must  be  here,  even  if  blank. 
[data  sets] 

6.2.  Control  Decks 

The  following  are  decks  of  UNIVAC  job-control  statements  each  identified  by  a  specific 
name.  These  control  decks  perform  various  functions  primarily  related  to  the  formating 
of  data  for  various  final  output  modes  (plots,  tables,  etc.).  From  the  user's  point  of 
view,  they  differ  from  the  programs  described  above  primarily  in  the  way  they  are  called 
from  mass  storage  and  executed. 

DIGIT 

DIGIT  is  a  deck  of  UNIVAC  job-control  statements,  stored  as  mass  storage,  which: 

(1)  executes  the  digitizer  program; 

(2)  takes  the  output  of  the  digitizer  program  and  inputs  it  into  CLEAN; 

(3)  executes  CLEAN; 

(4)  punches  the  output  of  CLEAN  complete  with  headers  and  labels 
(Caution:  The  labels  are  always  "PHOTON  ENERGY  (MEV)"  and 
"CROSS  SECTION  (MB)"); 

(5)  leaves  the  output  in  an  element  called  DIGOUT  for  use  in  other 
programs. 

The  use  of  DIGIT  is  best  illustrated  by  a  skeletal  example: 

@ELT,IS  DIGIN 

[digitizer  control  cards] 

@ASG,LJ0  IN., 8C, [user's  tape  number] 

@ADD,P  PF. DIGIT 

The  data  in  the  element  DIGOUT  may  be  accessed  immediately  in  the  same  job  by 
the  statement: 

@ADD,P  DIGOUT 

If  it  is  desired  to  use  the  data  in  the  same  job  and  not  punch  the  output  of  CLEAN,  precede 
the  @ADD,P  PF. DIGIT  card  with: 

SSETC  1 
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Two  digitizer  control  cards  must  be  given  for  each  curve  recorded  on  the  tape.  The 
order  of  these  pairs  of  cards  must  correspond  to  the  order  of  the  data  on  the  tape.  The 
first  card  of  each  pair  is  a  Header  Card;  the  second  a  Plot  Parameter  Card.  The  Header 
Card  gives  the  acquisition  number  of  the  data  set,  the  reference  number,  the  nuclide,  the 
reaction,  and  any  pertinent  comments  in  the  following  format: 


Position 


1 

-  5 

6 

7 

-11 

12 

13 

-18 

19 

-23 

24 

-30 

31 
32  -43 

44 
45 


52  -6 


Quantity 

Acquisition  number,  right  adjust 

Blank 

Number  of  data  points  (to  be  provided 

by  piogram) 
Blank 

Data  Center  reference  number 
Blank 

Nuclide/Element  identification 
25,25    Z,  right  adjust 
26,27    Symbol,  left  adjust 
28-30    A,  left  adjust 
Blank 

Reaction  with  (    ) ,  left  adjust 
Blank 
Source  flap;:  blank  if  digitized  at  Data 

Center;  or,  1  if  digital  data  was  supplied 

by  original  source 
Available  for  comments  if  needed 


The  Plot  Parameter  Card  must  contain  nine  numbers  separated  by  commas.  The  first  six, 
in  floating-point  notation,  give  the  coordinates  of  the  three  fiducial  points  used  to 
establish  the  axis  scales  for  the  digitizer  programs.  The  last  three  numbers  in  fixed-point 
notation  give  the  acquisition  number  for  the  plot,  an  x-axis  flag,  and  a  y-axis  flag.  For 
a  linear  scale  the  flag  is  0;  setting  the  flag  equal  to  1  indicates  that  the  corresponding 
axis  scale  is  logarithmic . 

Example: 


h- 
O 
Ld 
CO 

CO 
CO 

o 
cr 
o 


200« 

>. 

0  , 

r     i       i 

f 

1          V 
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Plot  Parameter  Card  entries: 


35 
ENERGY,     MeV 

15. ,0. ,35. ,0. ,15. ,200. ,932,0,0, 
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LIST 

LIST  is  a  combination  of  job-control  statements  (in  PP. LIST)  and  a  simple  program 
(in  PP.LISTR)  which  lists  data  in  the  breakpoint  file  (BF.)  on  the  slow  printer.  It  is 
designed  to  produce  camera-ready  copy.  Each  line  is  printed  twice  to  ensure  dark  type. 

Whenever  a  card  with  "$PAGE"  (starting  in  column  one)  is  encountered,  the  program  skips 
to  the  top  of  the  next  page  in  the  printer.  The  $PAGE  card  is,  of  course,  not  printed. 

Lines  longer  than  80  characters  may  be  printed  with  this  program.  The  first  80  char- 
acters of  the  long  line  are  entered  on  an  80-column  card  in  the  normal  fashion.  The  rest 
of  the  line  is  entered  in  columns  7-46  of  a  continuation  card  immediately  following  the 
primary  card.  The  continuation  card  must  have  "$IVDRE"  punched  in  columns  1-5.  A  maximum 
of  120  characters  may  be  printed  on  one  line.  The  primary  use  of  LIST  is  to  build  tables 
in  formats  suitable  for  mounting  on  H0%  reduction  paper.  The  long  linelength  required  by 
such  forms  precludes  the  use  of  simple  80-80  listing. 

There  are  two  uses  for  LIST: 

(1)  to  list  a  deck  of  cards: 

@ASG,T  BF. 
@DATA,I  BF.,,LISTR 
[card  deck  to  be  listed] 
@END  LISTR 
@ADD,P  PF.LIST 

(2)  to  list  an  already-existing  breakpoint  file  called  BF: 

@ADD,P  PF.LIST 
Example:  LIST  may  be  used  to  immediately  print  the  output  of  YTFLD4: 

@ASG,T  BF. 

@MAP,INX  MAPR, YIELDS 

IN  PF.  YIELDS 
@BRKPT  PUNCH$/BF 
@XQT  YIELD4 
[parameter  card] 
[data  sets] 
6BRKPT  PUNCH$ 
@ADD,P  PF.LIST 

PPLOT 

PPLOT  is  a  production  plotting  job-control  deck.  It  will  plot  points,  smooth  curves, 
or  zigzag  lines.  As  written,  the  program  produces  point  plots  in  a  standard  format. 
Normally,  the  nuclide  symbol  is  drawn  above  the  plot.  Also,  whenever  a  data  set  is  plotted 
in  POINT  mode,  the  reaction  and  reference  are  drawn  in  the  upper  right-hand  corner  of  the 
plot.  Two  control  cards  are  provided  for  altering  the  legend.  Provisions  are  made, 
however,  to  modify  this  format  both  through  the  use  of  plot  control  cards  as  well  as  by 
entering  the  desired  changes  for  specific  plot  parameters  which  control  trie   format  used  for 
a  given  plot.  An  arbitrary  number  of  data  sets  may  be  plotted.  Plotting  may  be  done  on 
the  printer,  the  SC  4020,  or  the  Calcomp.  It  is  recommended  that  a  printer  plot  be  made  on 
the  slow  printer. 
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A.  Program  Calling  and  Post-Processors 

The  plot  program  is  called  by  the  statement: 

@ADD,P  PF.PPLOT 

A  printer  plot  may  be  created  on  the  slow  printer  by  also  calling 

@ADD,P  PF.PRINTR 

The  normal  sequence  of  cards  is: 

@ADD,P  PP.PPLOT 

[control  cards] 

$$PARMS 

$  Parameter  Specifications  ^    (only  if  needed) 

$$END 

[data  sets] 

@ADD,P  PF.PRINTR  (if  printer  plot  is  desired) 

@MSG,W  [user's  account  number],  WRITE  ENABLE  REEL  [user's  reel  number] 

@ASG,TMJ  9.,8C,[user's  reel  number ]W 

gCLOSE   9 • 

@XQT  NBS*CALCMP. 

B.  Plot  Control  Cards 

PPLOT  is  controlled  by  a  set  of  control  cards.  On  all  of  these  cards,  a  dollar  sign 
must  appear  in  column  one.  The  control  cards  include: 


General  Controls 


$CURVE 


$LINES 

$POINT 
$CPLOT 


The  program  normally  plots  data  points.  This 
card  switches  to  a  curve-tracing  mode. 

This  card  switches  to  a  zigzag  line  tracing  mode. 
(This  is  also  useful  for  drawing  histogram  bar  graphs.) 

This  card  returns  the  program  to  point  plotting  mode. 

The  program  will  stay  in  CPLOT  mode  until  a  $FRAME 
card  is  encountered.  While  in  CPLOT  mode: 

-plots  are  12"  x  9.5" 

-the  nuclide  ID  is  suppressed 

-each  data  set  is  plotted  on  a  new  frame 

-all  these  frames  are  drawn  to  the  same  scale 

-the  data  sets  are  plotted  as  points  connected 
by  zigzag  lines 
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$SMALL         A  special  control  card  to  plot  small  frames.  This 
card  causes  a  4  x  2.5  in.  frame  to  be  plotted. 
Two  such  small  frames  will  fit  side-by-side  on  40% 
reduction  paper.  Or,  a  small  frame  may  be  pasted 
onto  a  regular  frame  as  an  insert.  The  effect  of 
$SMALL  persists  until  a  $FRAME  card  is  encountered. 

$FRAME         This  card  marks  the  beginning  of  a  new  frame.  The" 
program  returns  to  POINT  mode. 


Legend  Controls 

$LABEL         Anything  punched  in  card  columns  7  through  24  of 

this  card  is  written  as  a  single  line  into  the  plot 
legend. 

$PPLOT         This  card  causes  the  current  frame  to  be  repeated 

(re-drawn)  without  any  annotation.  The  primary  use 
of  this  card  is  to  plot  graphs  to  be  sent  to  a  drafts- 
man for  lettering.  If  entered  as  the  first  card  of 
a  frame,  $PPLOT  will  cause  a  "stripped"  frame  to  be 
plotted.  If  entered  as  the  last  control  card  of  a 
frame,  two  plots  will  be  drawn,  one  with  annotation 
and  one  without. 

Data  Modification  Controls 

The  energies  and/or  cross  section  values  of  data  sets  may  be  renormalized  by  the 
control  cards: 

$MUL  E[value]       Multiplies  the  energies  of  the  next  data  set  by  the 

given  floating-point  value. 

$MUL  S[value]       Multiplies  the  cross  sections  of  the  next  data  set 
by  the  given  value. 

Similarly,  offsets  can  be  added  to  the  data  by: 

$ADD  E[value]       Adds  the  value  to  the  energies  of  the  next  data  set. 

$ADD  S[value]       Adds  the  value  to  the  cross  sections -of  the  next- data 

set . 

All  four  of  the  above  manipulate  only  the  next  data  set.  If  used  in  combination,  the 
multiplication  is  done  first.  All  values  are  given  as  floating-point  numbers  which  must 
start  immediately  after  the  E  or  S  (i.e.,  no  space). 

C.  Parameters  and  Defaults 

PPLOT  does  not  require  any  plot  parameters  to  be  specified.  The  normal  plot  constants 
are: 

XDIST  =1.25 
YDIST  =1.25 
XLNGTH  =9.25 
YLNGTH  =6.00 

FONTB  =  .09375  -FONT  is  plot  symbol  size.   .09375  =  3/32  in. 
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YDOWN 

SYMBOL 

YUP 

XSTART 

YDIV 

YSTART 

YDECPL 

FONTH  =  .09375 

CHARB  =  0.1875  -CHAR  is  legend  and  numeric  symbol  size. 

CHARH  =  0.1875  0.1875  =  3/16  in. 

SCALB  =  0.15625  -SCAL  is  axis  label  symbol  size. 

SCALH  =  0.15625  0.15625  =  5/32  in. 

XLBL  =  PHOTON  ENERGY  (MEV)      -If  only  one  data  set  is  plotted,  XLBL  and 

YLBL  =  CROSS  SECTION  (MB)        YLBL   are  coPled  f*rom  the  label  cards  of  the 

data  set. 

The  plot  variables  are  normally  calculated  from  the  data.  They  include 

XDOWN 
XUP 
XDIV 
XDECPL 

The  various  plot  parameters  and  variables  are  defined  in  the  figure  given  on  p.  21. 

The  plot  control  switches,  and  defaults,  are: 

RULE  =  1  Set  to  zero  to  inhibit  upper  right  legend. 

NTOP  =  1  Set  to  zero  to  inhibit  nuclide  ID. 

MOVE  =  1  Set  to  zero  to  inhibit  frame  advance. 

RECALC  =  0  Set  to  1  to  re-calculate  XDIV,  YDIV,  XDECPL, 

YDECPL,  XSTART,  and  YSTART. 

Any  and  all  of  these  parameters  (in  fact,  the  entire  SPECS  array)  can  be  modified  by  means 
of  the  $$PARMS  cards.  The  constants  and  variables  are  meant  to  be  changed.  However,  it  is 
recommended  that  the  casual  user  leave  the  plot  switches  alone. 

Parameters  are  entered  in  FORTRAN  "NAMELIST"  format.  That  is,  you  punch  the  parameter 
name,  an  equals  sign,  the  numeric  value,  and  a  comma. 

Example : 

$$PARMS 

$    XDOWN=10 . ,XUP=40 . ,XDTV=4 . ,XDECPL=0 . , 

$$END 

Don't  forget  the  dollar  signs  in  column  one.'  All  parameters  return  to  their  default  values 
when  "$FRAME"  is  encountered. 

D.  Examples : 

Actual  use  of  PPLOT  is  best  illustrated  by  a  few  examples. 

1.  To  plot  a  (y,sn),  (y,2n),  and  (y,p)  graph  in  a  form  ready  for  the  draftsman, 


use: 


@ADD,P     PF. PPLOT 
[(y,sn)  data] 
C(Y,2n)  data] 
[(Y,p)  data] 
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$CURVE 

[(Y,sn)  curve] 

[(Yj2n)  curve] 

[(y,p)  curve] 

$PPLOT 

@ADD,P  PF.PRINTR 

@MSG,W  [user's  account  number],  WRITE  ENABLE  REEL  [user's  reel  number] 

@ASG,TMJ  9 .,8C, [user's  reel  number ]W 

@CLOSE  9 . 

@XQT  NBS*CALCMP. 

This  will  give  both  a  printer  plot  and  a  tape  for  the  Calcomp.  To  do  the  same  thing,  but 
with  a  nine-inch,  rather  than  a  six-inch  plot  height,  use: 

@ADD,P  PF.PPLOT 

$$PARMS 

$    YLNGTH=9 . , 

$$END 

[(YjSn),  (v,2n),  and  (y»p)  data  sets] 

$CURVE 

[(y,sn),  (y,2n),  and  (y,p)  curves] 

$PPLOT 

@ADD,P  PF.PRINTR 

etc. 

2.  To  plot  a  (Y»xn)  cross  section  and  a  (y5p)  cross  section  on  separate  frames,  use: 

@ADD,P  PF.PPLOT 

[(y,xn)  points] 

$CURVE 

[(y,xn)  curve] 

$FRAME 

[(y,p)  points] 

$CURVE 

[(y,p)  curve] 

@ADD,P  PF.PRINTR 

etc. 

3.  To  plot  a  (y»p)  cross  section  and  a  (Y>Mt)  cross  section  on  separate  frames,  with 
the  intention  of  pasting  the  (y,u.)   graph  in  the  upper  right-hand  corner  of  the  (y>p)  graph 
as  an  insert,  and  also  assuming  tne  (yjM^)  data  consist  of  bar  graph  corners,  use: 

@ADD,P  PF.PPLOT 

$$PARMS 

$        XSTART=1.5,  (this  moves  the  legend  to  the  upper  left  comer) 

$$END 
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[(y,p)  points] 

$CURVE 

[(y,p)  curve] 

$PPLOT 

$FRAME 

$SMALL 

$LABEL  (G,MU-T)  64ER2 

$LINES 

[(Y,ut)  data] 

$PPLOT 

@ADD,P  PF.PRINTR 

etc. 

4.  To  obtain  a  set  of  comparison  plots,  use: 

@ADD3P  PF.PPLOT 

$CPLOT 

[all  (y,n),  (y,sn),   (y,xn),  (y,vt)   data] 

$FRAME 

$CPLOT 

[all  (y>p)  and  (y5xp)  data] 

$FRAME 

$CPLOT 

[all  (Y,yt)  data] 

$FRAME 

$CPLOT 

[all  (Yjc)  data] 

etc. 

SPCDST 

This  deck  creates  an  index  table  of  spectrum  and  angular  distribution  papers.  Use  is 
made  of  the  Photonuclear  Data  Index  which  must  be  in  mass  storage  in  a  file  called 
"DC*INDEX.".  The  output  table  is  placed  in  a  file  called  "BF.",  where  it  can  be  listed  by 
LIST.  The  user  enters: 

1)  the  nuclide 

2)  the  references  and  reactions  of  SPC  entries 

3)  DST 

4)  the  references  and  reactions  of  DST  entries 

5)  @EOF 

The  references  must  start  in  Column  1.  The  user  should  leave  at  least  one  space  after  the 
reference  as  a  separator.  The  reactions  do  not  need  parentheses. 
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Calling  sequences  and  example: 

@ASG,T  BF. 

@ADD,P  PF.SPCDST 

5B  10 

68SH6  G,P 

DST 

68SH6  G,P 

63FU3  P,G 

65PA1  HE,G 

68LI1  HE,G 

@EOF 

@ADD,P  PF.LIST 
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